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Search for SUSY at LHC in the first year of 

data-taking. 

Abstract 

If Supersymmetry would manifest itself at a low mass scale it might be 
found already in the early phase of the LHC operation. Generic signatures for 
Supersymmetry in pp-collisions consist of high jet multiplicity, large missing 
transverse energy as well as leptons in the final state. The CMS search strategy 
and prospects for a SUSY discovery in the first year of data-taking is reviewed. 
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1 Introduction 



The Large Hadron Collider (LHC) [lj at CERN opens a new energy regime offering 
a very exciting discovery potential for physics beyond the Standard Model. In this 
paper the search for Supersymmetry (SUSY) in the startup scenario at yfs = 10 TeV is 
discussed. Supersymmetry exists in many theories beyond the Standard Model. Two 
t-h of the best studied braking mechanisms are the Minimal SuperGravity (mSUGRA) 

Psj and Gauge Mediated Symmetry Breaking (GMSB). Benchmark points have been 

denned in the framework of mSUGRA and GMSB to study various experimental 
SUSY signatures. In the CMS experiment [2] SUSY analyses are organized according 
to topologies, e.g. number of leptons, photons and jets in the final state, which arise 
from SUSY cascades. A detailed definition of the benchmark points can be found 
here j3]. 

The paper is organized as follows: in Section [2] the SUSY search in multijet final 
state is discussed [3]. Searches for new physics in diphoton jS] and dilepton [6] final 
state are described in section [3] and [3] respectively. 

2 Search for SUSY in fully hadronic final state 

This section describes a search strategy for a possible discovery of SUSY signatures 
with the CMS detector at the LHC using exclusive n-jet events (n = 2. . .6). The 
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event topology under investigation consists of n high-px jets and two invisible neu- 
tralinos which lead to a missing energy signature. The high-pp jets are produced in 
the decay chains of the initially produced heavy squarks and gluinos. The main aim 
of the analysis is to develop a robust measurement technique suitable for the early 
physics data at the LHC and stable with respect to jet energy mismeasurements. 
Before applying any event selection, multijet production from QCD is the dominant 
process, where missing energy is introduced through jet mismeasurements. 

2.1 Trigger and selection 

The benchmark points of mSUGRA LMO and LM1 are used to estimate the trigger 
efficiency for signal events. Both signal points have 100% efficiency after all cuts for 
the single jet trigger HLT_JetllO (one jet with corrected jet transverse momentum > 
110 GeV/c). 

Hadronic jets are reconstructed from calorimeter energy deposits which are clus- 
tered using an iterative cone algorithm with R = 0.5 [7]. Furthermore, these jets 
are required to have a transverse momentum greater than 50 GeV/c, pseudo-rapidity 
1^1 < 3.0, and an electromagnetic fraction F em < 0.9. The transverse momentum 
of the leading jet and second leading jet need to exceed 100 GeV/c and the pseudo- 
rapidity of the leading jet is required to be smaller than two. Based on the jets defined 
above two additional variables are defined: Ht as the scalar sum over the transverse 
momenta of the selected jets in an event, Ht = YliPr an d the missing transverse mo- 
mentum of the event calculated as H™ ss = — X^Pt'' 1 - In order to reduce background 
events from SM processes Ht is required to be greater than 350 GeV/c. 
All events are rejected where either an isolated muon [8] with px > 10 GeV/c or an 
isolated electron [U] with p T > 10 GeV/c or photons [3] with p T > 25 GeV/c or jets 
with p T > 50 GeV/c that does not fulfill the other criteria (\r}\ < 3 or F cm < 0.9) are 
found. 

2.2 Analysis method and results 

In the following a kinematic variable («t) is used that allows separation of signal 
events with real missing energy from QCD events in which missing energy is created 
by jet energy mismeasurements. 

In the dijet case (n = 2) transverse momentum conservation requires the px of the two 
jets in QCD events to be of equal magnitude and back-to-back in the azimuthal angle 
cf). The variable op, first introduced in Ref. [TO], exploits exactly this requirement. It 
is defined as 

a T = E J T 2 /M T (1) 
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Table 1: Numbers of events expected for the dijet case and for n = 3 . . . 6 jets case 
for background samples (QCD, Z — > vv+jets, W+jets, tt and Z+jets) and the LMO 
and LM1 signal points. The final numbers of selected events are shown after the cuts 
on a T and R(H™ SS ). 



where E^ is the transverse energy of the second leading jet in the event and Mt is 
defined as 
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and n = 2 in the dijet case. For a well measured QCD dijet event, E^ = 0.5 x 
and H™ 1SS = 0, thus is exactly 0.5. 

To define for more than two jets the n-jet system is reduced down to a two- 
jet system by combining jets into two pseudo-jets. The E T of the pseudo-jets is 
calculated as the scalar sum of the contributing jet Et- All possibilities of how n jets 
can be combined into two are tested and the combination is chosen where the resulting 
pseudo-jet are most similar, i.e., for which the difference Aifr = E™ 1 — E™ 2 is 
minimal. For n jets, «t is then obtained in the same way as in Eq. [T] 

In the event selection is required to be greater than 350 GeV/c which is well 
above the transverse momentum threshold of 50 GeV/c for a single jet. However 
several jets below that threshold could still lead to a considerable amount of ignored 
momentum in the event. For that reason the if ™ 1SS determined using all jets having 
a pt larger than 30 GeV/c , if™ lss (jetpr > 30 GeV/c), is calculated and compared to 
the if™ lss determined from the selected jets only, if™ lss (selected jets). The ratio 

R(H™ iss ) = H™ iss (selected jets) /H™ iss (jet p T > 30 GeV/c) (3) 

can be used to single out events where the inclusion of lower momentum jets does 
significantly improve the balance of the event. If the missing transverse energy (if™ lss ) 
is increased by 25% due to the fact that the transverse momentum threshold of the 
selected jets is 50 GeV/c and not 30 GeV/c the event is rejected, thus R(H™ ms ) is 
required to be smaller than 1.25. 

The «t distributions for the dijet case and the sum of n — 3 ... 6 jets case are 
shown in Fig. [l] where the requirement on R(H™ ms ) has already been applied. In 
both figures the QCD background peaks, as expected, sharply at a value of 0.5. To 
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(a) Distribution of ay for dijet events. (b) Distribution of ay for events with n = 

3 ... 6 jets. 

Figure 1: «t distribution. 

account for finite jet energy and resolutions events are only selected if «t is larger 
than 0.55. 

The resulting event yields for signal and background are summarized in Tab. [TJ 
All expected event yields correspond to an integrated luminosity of 100pb _1 . It can 
be seen that in the dijet case only Z — > vv + jets and W + jets events give a small 
background contribution over a clear signal. At higher jet multiplicities n — 3 . . . 6, 
top decays as well as about one QCD event contribute to the remaining background 
after the final selection. 



2.3 Establishing a Signal incompatible with Standard Model 
Background in Data 

To establish the discovery of a SUSY signal the fact that signal events are produced 
more centrally in pseudo-rapidity compared to the SM backgrounds, in particular 
compared to QCD events whose main production mechanism is t-channel exchange, 
is used. The pseudo-rapidity of the leading jet can be used as a measure of the 
centrality of an event. For the SM background the ratio i?„ T (0.55) of events with 
larger than the cut value over that of events with «x smaller than the cut value, is, 



as shown in Fig. 2(a), approximately constant as a function of pseudo-rapidity and 



independent of i/j- -R aT (0.55) behaves very differently in the presence of a SUSY 



signal as illustrated in Fig. 2(b) for the LMO benchmark point. As can be seen, 
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(a) SM background only 
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(b) SUSY(LMO) + SM background 



Figure 2: i? aT (0.55) as a function of \r]\ for different cuts for SM background 
only and in case of presence of a SUSY signal. The error bars indicate the expected 
statistical uncertainties for a data sample of 100 pb _1 . 



the presence of a SUSY signal manifests itself with two distinct features: R aT (0.55) 
exhibits a negative slope with larger values of \rj\; tighter requirements on Ht result 
in a steeper slope and an offset in i? QT (0.55). For events with Ht > 350 GeV/c the 
measured i? QT (0.55) in the central \r)\ bins is well above the ratios obtained from the 
control region 300 < < 350 GeV/c and increases with smaller values of 77. Even 
with a systematic uncertainty of 100% on R aT (0.55) in the control region the excess 
would remain convincing. 



3 Search for SUSY in diphoton final state 

The final state with two high Et photons and large missing transverse energy can 
harbor new physics signals in a variety of theoretical scenarios, most notably GMSB. 
In this section, a data-driven strategy that can be used in a plausible start-up scenario 
to predict the E™ 1SS distribution in a diphoton sample from the SM processes is 
described. Observation of an excess of events at high E™ 1SS would be a signature of 
new physics. 

The SM contribution to diphoton plus E<f lss final state is small. The only physics 
backgrounds are Z77 — > z/z/77 and W77 — > £^77. The instrumental background 
has three major components. The first component results from QCD events with no 
true E™ 1SS (QCD background), such as direct diphoton, photon plus jets, and multijet 
production. The second component comes from events with real E™ ms (Electroweak 
background). This component is dominated by Wj and Wj production where the W 
decays into an electron plus a neutrino when the electron is mis-reconstructed as a 
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photon. The third background is associated with the high energy muons from cosmic 
rays or beam halo (Non-beam background). 



3.1 Trigger and selection 

The trigger applied in this analysis requires a photon with E^ >25 GeV. This trigger 
is expected to run unprescaled for the duration of the first CMS run, and is fully 
efficient for the SUSY signal considered. 

For the photon candidate selection, two objects in the ECAL Barrel (|?7| < 1.45) 
with pt >30 GeV/c were required to be isolated and to have a negligible deposit 
of energy in the hadronic calorimeter around the photon. The selected objects are 
classified as electrons if they have an associated track stub in the pixel detector, 
referred to as pixel seed, and as photons otherwise. In this way three independent 
samples are defined: 77 sample comprising events with at east two photons; ej 
sample with at least one electron and at least one photon; ee sample with at least 
two electrons. 

Finally a cut E™ ss > 80 GeV is applied to reduce significantly the SM background. 
The results of the selection are reported in Tab. [2] 





77 


7+jet 


multijets 


W, W+n-f 


Z, Z+n^f 


GMlc 


all events 


1055 ± 4 


3189 ± 100 


173 ± 37 


8.5 ± 3.0 


23± 1.3 


20.7 


after E™ iss 


1.3 ± 0.16 


1.3 ± 0.16 




0.09± 0.04 


0.09± 0.02 


14.8 



Table 2: Event counts before and after the cut on E^^ for the 77 sample for 100 pb 
at 10 TeV. 



3.2 Non-beam background 

Energetic cosmic muons or muons from beam halo can emit photons as they pass 
through the ECAL. The shower shape and arrival time of these photons is slightly 
different from the ones originating from the interaction point, especially for beam 
halo. 

Photons from the interaction point hit a crystal approximately perpendicular to 
the face because the crystals are rotated in rj such that they point back to the in- 
teraction point. As a result of this rotation, photons from the interaction point give 
rise to showers with a small spread in 77. However, photons from the beam halo 
traveling parallel to the beam tend to have a large spread in r]. Cosmic ray photons 
come from all possible angles and could have very narrow or very wide showers. Fig. 
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Figure 3: Photon shower and Timing in beam halo, cosmics, and photon gun events. 



3(a) displays the distribution of the energy- weighted RMS of the shower profile in rj 
(77-width) for the three sources of photons. 

The expected time resolution for individual crystals is better then 1 ns, and the 
non-beam backgrounds can be discriminated against using the measured cluster time. 
Cosmic background is asynchronous, but beam halo has a very specific time distri- 
bution. Protons and halo muons travel parallel to each other and have the same 
z(t). The time difference between the photons from the IP and halo is then given 
by At = (Z + \J Z 2 + E?)/c. The value of R, the radius of the shower in the ECAL, 
is not known. Fortunately, the resulting uncertainty in At is fairly small, as shown 



in Fig. 3(b) , which shows the comparison of the measured time of the photons from 
beam halo, cosmics, and prompt photons, corrected for the detector geometry. 

Although the exact amount of non-beam backgrounds is very hard to predict, a 
combination of the two techniques explored above (shower shape and timing) will be 
sufficient for the effective elimination of these backgrounds. It will also be possible 
to isolate a pure sample of non-beam backgrounds, ensuring a reliable estimation of 
the residual contamination. 



3.3 Determination of E™ ss Distribution from Data 

The contribution of the Electroweak background is determined using the ej sample. 
The E™ 1SS distribution in sample is multiplied by / e _» 7 /(l — /e-»- 7 ), where / e _^ 7 
is the electron-photon mis-identification rate. /e-> 7 is obtained by fitting the mass 
of the Z for the ee sample and e7 sample, and by comparing the integrals of these 
fits. The resultant distribution is taken as an estimate of the this background and is 
subtracted prior to assessing the background from QCD events. 
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(b) SUSY(GMlc) + SM background 



Figure 4: Background closure test using Z — > e + e events to describe the QCD 
background. 



The E™ 1SS distribution for the QCD background is modelled using a sample of Z— > ee 
events, selected from ee sample with invariant mass cut on the two electrons, 80 
< M ee < 100 GeV. 

The key assumption is that the di-EM system which is measured comparitively 
well, and the recoil and other hadronic activity which is measured poorly can be 
separated in the event. The di-EM is used as a measure of this additional activity. 
Its amount is different for the Z— ► ee and the 77 events. To obtain the proper shape 
of the E™ 1SS distribution, the Z— > ee events are reweighted so that their di-EM px 
distribution matches the one in the 77 sample. The weighted E™ 1SS distribution is 
then normalized to the 77 sample in the low missing transverse energy (E™ 1SS < 
20 GeV) region. 

Fig. 4(a) show the closure test for Z — > e + e~ sample used for QCD background 
determination in absence of SUSY signal. The presence of the signal can bias the 
background estimation, thus the entire analysis is repeated with GMlc SUSY signal 



mixed in, as shown in Fig. 4(b) 



The data driven estimation of the background agrees very well with the number 
of expected events evaluated with the MC truth, as shown in Tab. [3] 



4 SUSY discovery potential and measurement of 
a dilepton mass edge 

In mSUGRA models the lightest neutralino escapes detection and no mass peaks can 
be observed in SUSY decay chains. Of special interest are robust signatures such as 
edges in mass distributions in leptonic final states which can be probed with the CMS 
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ivrQCD 


atEWK 
iV 77 


jyGMlc 




N EWK 


no SUSY 


2.61 ± 0.23 


0.17 ± 0.04 




2.34 ± 0.65 


0.35 ± 0.10 


with SUSY 


2.61 ± 0.23 


0.17 ± 0.04 


14.8 ± 0.1 


2.48 ± 0.67 


0.50 ± 0.10 



Table 3: Closure test for 77 sample with and without GMlc SUSY signal. The 
number of events corresponds to 100 pb" 1 at 10 TeV. 



experiment. 

The purpose of this analysis is to observe a significant excess of opposite sign 
same flavour leptons over the various backgrounds and to measure the endpoint in 
the invariant mass distribution. All flavour symmetric background (including SUSY 
decays of this type) can be determined from data events with opposite sign opposite 
flavour leptons. The aim is to perform such an analysis already with the first LHC 
data which is expected to amount to roug hly 200-300 pb" 1 in 2010. 

The leptonic decay of the next-to-lightest neutralino gives a characteristic signa- 
ture. This decay can proceed in different ways even in the mSUGRA model. A mass 
difference of the neutralinos smaller than the Z boson mass and any slepton mass 
leads to a three body decay. In that case the endpoint in the lepton invariant mass 
represents directly the mass difference of the two lightest neutralinos 

mu,max = m^o - m^o (4) 

A two-body decay occurs via a real slepton and is allowed if at least one slepton is 
lighter than the mass difference of the neutralinos. In that case the endpoint can be 
expressed by 

(m 2 - mf) (mf - m| ) 

\ m U ) - — 2 \°) 

where mj is the mass of the intermediate slepton. The shape of the mass edge results 
only from kinematics and is triangular. 



4.1 Trigger and selection 

Single lepton (electron or muon) trigger is used in this analysis. Since the leptons 
originating from the signal decay have a very soft pt spectrum, the triggers with the 
lowest available threshold for electrons (px > 15 GeV/c) and muons (pt > 11 GeV/c) 
are used. 

The base selection requires two isolated leptons of opposite sign. Muon identifica- 
tion requires reconstruction in both the muon system and the inner tracker [S] . Each 
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Figure 5: DiLepton invariant mass. The black solid line represents the extrapolation 
from the opposite flavour distribution. The black points represent pseudo data of one 
experiment where no scaling has been applied but exactly 200 pb _1 of MC events 
have been analysed. 



electron has to fulfill the tight electron identification criteria, which consist of a set 
of cuts depending on the electron px and rj j9]. Additionally a p-r > 10 GeV/c and 
1 77 1 < 2 is required for each lepton. To stay above the trigger threshold the first lepton 
is required to have a > 16 GeV/c. The main SUSY selection is based on jets and 
missing transverse energy. The cuts have not been optimized at a certain benchmark 
point, but should reflect the general SUSY signature. The selection requires three 
jets with > 100 GeV/c, p!£ > 50 GeV/c, and pj, > 50 GeV/c. A missing transverse 
energy of at least 100 GeV is required. 



Sample 


tt+jets 


Z+jets 


W+jets 


Diboson 


Dijet 


LMO signal 


LMO incl. 


■N events 


80 


1 


2 








87 


362 



Table 4: Number of selected events using the described event selection for an inte- 
grated luminosity of 200 pb -1 . 



The number of events obtained after the selection is listed in Tab. HI All back- 
ground which leads to uncorrelated lepton pairs can be measured directly from data [H] 
In order to extrapolate to the same flavour opposite sign lepton pair distribution the 
opposite sign opposite flavour lepton pairs and use this distribution are selected. With 
this method one is able to predict all backgrounds which produce uncorrelated leptons 
such as W, it, dijet and WW events. 

The invariant mass distribution of all opposite sign same flavour leptons for 
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200 pb 1 is shown in Fig. 5(a) The opposite sign opposite flavour distribution used 
to extrapolate the background is displayed in Fig. 5(b) 



4.2 Determination of the mass edge 

The model used for the fit of the mass edge consists of three parts. To model the 
signal the theoretical model [121 EE] convoluted numerically with a gaussian has been 
used in case of a 3-body decay 



^y A -y 2 (m 2 + M 2 ) + (m M) 2 



2rca J (y 2 — m|) 2 



o 



-( m u~y) 2 



x (-2y 4 - y 2 [m 2 + 2M 2 ) + (m M) ) — (6) 

where m = m^o — m^o is the difference, M = m^o + m^o is the sum in neutralino 

A2 Al 7 A2 Al 

mass and Mz is the Z mass, which is kept fixed. In case of the two-body decay the 
signal model consists of a triangle convoluted with a gaussian 

rricut 

T(m u ) = / dy ■ ye & (7) 

Vzna J 
o 

A curve parametrized as 

B(m u ) = ml ■ e- b - m " (8) 

has been used to fit the opposite sign opposite flavour invariant mass distribution. 
Additionally the Z peak is fitted using a Breit-Wigner convoluted with a gaussian. 
The number of signal N$i g , background Nskg and Z events Nz are fitted as well. 

4.3 Expected results 

A simultaneous fit to the the ee, /i/i (signal plus background model) and e/i (back- 
ground model) invariant mass distributions is performed. 



The fit to the invariant mass distributions at LMO is shown in Fig. 6(a) and 
it yields a value of mu ;max = (52.8 ± 1.1) GeV for the dataset of exactly 200 pb -1 
and mu^max = (51.3 ± 1.5) GeV for the full MC data where the error is the one 
expected after 200 pb -1 . The derived number of signal events n sig = 97 ± 17 agrees 
with the number of signal events from MC truth (Tab. [4]). The theoretical endpoint 
nT>n,theo = 52.7 GeV is reproduced in case of the fit with the three-body decay model. 
The background fit of the opposite sign opposite flavour lepton pairs is shown in 



Fig. 6(b) The total number of background events is n^kg = 201 ± 14 which is in 



agreement with the expected number from MC truth (192). 
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(a) ee, \x\i invariant mass distributions 
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Figure 6: The combined fit at LMO, the green curve represents the SUSY signal 
model, the red curve is the background function and the light green dashed line the 
Z contribution. 

The main sources of systematic uncertainties are the jet energy scale, the electron 
energy scale, the lepton efficiency and the modeling of the background and of the 
resolution. Combining systematic and statistical errors, the expected results for LMO 
and 200pb~ 1 at 10 TeV is 

m u ,max = (51.3 ± 1.5 stat . ± 0.9 syst .) GeV/c (9) 

compared to a theoretical value of 52.7 GeV. 



5 Summary 

The strategy for three different analyses aimed to discovery SUSY in the first year 
of data-taking at LHC considering a/s = 10 TeV with the CMS experiment were 
presented. 

The search for SUSY in multijet final state is carried out in the context of SUSY for 
several sets of parameters in the mSUGRA parameter space assuming an integrated 
luminosity of 100 pb" 1 . The discrimination power of «t against SM background from 
QCD events provides, for favourable SUSY benchmark points, signal over background 
ratios of 4 to 8 depending on the considered jet multiplicity bin. 

In the context of the GMSB models a strategy to infer the existence of new physics 
in diphoton events with large missing transverse energy has been developed. Back- 
ground contribution in the high E!f lss region is estimated can be precisely estimated 
with a pure data-driven method in 100 pb -1 . 
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A significant excess of SUSY opposite sign same flavour lepton pairs can be found 
within the first 200 ptT 1 at LMO. The signal provides a quite robust signature and 
the background determination directly from data is possible. At LMO the combined 
fit of the dileptonic endpoint is possible with 200 pb -1 with an expected uncertainty 
of 1.8 GeV. 
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